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HERMITE-HADAMARD-TYPE INEQUALITIES IN THE APPROXIMATE 

INTEGRATION 
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£N| ■ ABSTRACT. We give a slight extension of the Hermite-Hadamard inequality on simplices 

and we use it to establish error bounds of the operators connected with the approximate 
integration. 



1. Introduction 



in 

(N 

"^ , In the numerical analysis there are many classical quadrature rules. For each of them 

\J—{ ■ the error terms are well known. In the series of recent papers J5] |4j [6) the author estab- 

C"| , lished error bounds of quadrature operators under regularity assumptions weaker from the 

classical ones. The method was based on the use of Hermite-Hadamard-type inequalities 

for convex functions of higher order. 

In this paper we develop an analogous idea for convex functions of several variables 
defined on simplices in R". It is possible because the generalizations of the Hermite- 
Hadamard inequality are known for multiple integrals, e.g. by Choquet theory. Recently 
Bessenyei proved in [1J such an inequality on simplices using an elementary approach. 
2^ ' We slightly extend this result and we use this extension to give error bounds of operators 

connected with the approximate computation of multiple integrals (i.e. the cubatures). 

Let po, . . . ,p n e R" be affine independent and let S be a simplex with vertices po, . . . , p n . 
Denote by p its barycenter, i.e. 

_ 1 A 

i=0 

and by vol(S) its volume. 

In our proofs we use some results contained in [ 1 1 . One of them is the following 



Lemma 1. If A : R™ — ► R is an affine function then 

A{p) = —^ f A{x)dx = -}—JTA{ Pi ). 
vol(S) J s n + 1 7^ 

The properties of convex functions needed in this paper can be found in many textbooks 
and monographs, for instance in [2| (cf. also HI). 

Recall that a linear functional T defined on a linear space X of (not necessarily all) 
functions mapping some nonempty set into R is called positive if 

/ < g => T(f) < T(g) 
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for any /, g E X. 

2. Hermite-Hadamard-type inequalities 
We start with a slight generalization of the Hermite-Hadamard inequality. 

Theorem 2. Let T be positive linear functional defined (at least) on a linear subspace of 
all functions mapping S into K generated by a cone of convex functions. Assume that 

vol(S) J s 

for any affine function A : S — > R. Then the inequality 

n 

n+l t J o 
holds for any convex function f : S — ► R 

Proof. By convexity the subdifferential df(p) is non-empty, whence there exists an affine 
function A : S — > R supporting / at the point p (i.e. A(p) — f(p) and A ^ / on 5). 
Using Lemma[T]and the properties of T we get 

/(F) = A(p) = -^-r f A{x)dx = T(A) < T(/). 

V01(£>J Js 

Since po, ■ ■ ■ , P™ are affine independent, there exists (exactly one) affine function B : S — > 
R such that B(j>i) — f(pi), i — 0, . . . ,n. A simplex 5 is the convex hull of the set 
{po, . . . ,p n }- Then by convexity / ^ B on S. Using once again Lemma Q] and the 
properties of T we arrive at 

nf)^nB) = ^J s B( X ) dx = ^ l ±Bi Pi ) = ^ l tf( Pi ). 

a 

As an immediate consequence we obtain the Hermite-Hadamard inequality. 
Corollary 3. Iff:S^>Ris convex then 



1 /• 1 " 

vol(S) 7 S n + 1 j-^ 



Another important consequence of Theorem [2] concerns the operators of the form of a 
convex combination of the values of / at some appropriately chosen points of S. 

Corollary 4. Let m € N, Ai, . . . , A TO ^ wirt X)I=i A* = 1 ant/ £ 1( . . . ,£ TO e S 1 vvj'fn 
Sj=i ^»£« = P- If a function f : S -tKis convex f/zen 



m -. n 

/(P)<£a,/&K— --^/(k). 



i=0 

Proof. Clearly 



?"(/) :=!>/(& 



1=1 
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is a positive linear functional. Let A : S — > R be an affine function. Then by LemmaQ] 

1 



t(a) = e A ^fe) = MJ2 A ^J = A (p) 



vol(S) 



A(x)dx 



and all the assumptions of Theorem |2] are fulfilled. 



□ 



This result is useful in the approximate computation of multiple integrals over S. To- 
gether with further results of this paper it allows us to give error bounds of the cubatures of 
the form of a conical combination of the values of an integrand at some points of S. These 
error bounds depend on the second order differential of the function involved. 

To conclude this section notice that in Theorem[2]not a form of an operator T is essential 
but three its properties are important: T is linear, positive and it coincides with the integral 
mean value in the class of affine functions. 



3. A NORM OF A SECOND ORDER DIFFERENTIAL 

Let ip : R™ — + R be a quadratic form, i.e. 

n 

<p(x) = E a ij x l Xj 

for a: = (x\, . . . ,x n ) S R n . Let 

IMI :=sup{|(^(a;)| : ||a;|| = l}. 

This is well defined norm on a linear space of all quadratic forms of n variables with the 
following properties: 



MK E 



H] | j 



*iJ=l 



|^)| < (Ml -N 2 , xel" 

The details may be easily checked (cf. e.g. Q). 

Now let / e C 2 (S). Then the partial derivatives of the second order are bounded, say, 
there exists a constant M > such that 

d 2 f 



dxidxj 



(«) 



< 



M 



for any u e S, i,j = 1, . . . ,n. Then for any u G S 

1,3 = 1 J 



(u)ViVj 



whence 



This allows us to define 



l<*7(«)ll < E 



o 2 f 



dxidxj 



(u) 



< M. 



K/Hoo := sup d 2 /(u) 

m6S 
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Obviously 

(1) 
(2) 
(3) 



\\d 2 (af)\\ oo = \a\-\\d 2 f\U 
K/(u)(t; ) t;)|<||d 2 /(«)||.|H| 2 <||d 2 /|UH| a 



for any f,g€ C 2 {S), a G R, u G S and v G R™. 



Lemma 5. /// G C 2 (S) and g(x) 
convex. Moreover, 



_ IM 2 /llc 



f/ien the functions g + f and g — f are 



\\d 2 (g + f)\\ oo ^2\\d 2 f\U and \\d 2 (g - f)]^ ^2\\d 2 f\\ oc . 

Proof. It is easy to compute d 2 g(u)(v, v) — ||d 2 /||oo II'L'H 2 - Then by (0 

\d 2 f(u)(v,v)\ ^ ||d 2 /l|oo||«|| 2 = d 2 g(u)(v,v). 

Hence, for any u G S, v G R™, 

d 2 (g + f)(u)(v,v)^0 and d 2 (g - f)(u)(v,v) > 0, 

which means that for any u G S , d 2 (g + f ) (u) and d 2 (g — /) (u) are positive semi-definite. 
Then g + / and g — / are convex. 

If ||t;|| = 1 then d 2 g(u)(v,v) = ||rf 2 /||oo f° r an Y « G S 1 , which implies ||d 2 g(u)|| = 
\\d 2 f\\oo,u G 5. Then ||<i 2 g||oo — ||gP/||oo and by (Q}, (O we obtain the second assertion 
of the Lemma. □ 



4. Applications to the approximate integration 
Lemma 6. Iff G C 2 (S) then 



f(x)dx-vo\(S)f(p) 



^IM 2 /Hoo f \\x-p\\ 2 dx. 



Proof. Write Taylor's formula with the remainder depending on the second order differen- 
tial: for any ieS there exists a point £ x belonging to the segment with endpoints x and p 
such that 



(4) 



f(x) - m - df(p)(x-p) = ±d 2 f(z x )(x~p, x -p). 



Since the left-hand side is an integrable function of x G S, then so is the right-hand side. 
Observe that x i— ► df(p)(x — p) is an affine function, whence by LemmaQ] 



df(p)(x - p)dx = df(p)(p - p) = 0. 



Therefore integrating both sides of © and using (0 we get 



f(x)dx - vol(S)/(p) 



d 2 f(£ x )(x - p, x - p)dx 



1 
^2 



\d 2 f(^)(x-p,x-p)\dx^-\\d 2 f\\ c 



\x — p\\ dx. 



a 
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Theorem 7. Let T be a positive linear functional defined on C 2 (S) such that T(p) 
fa p(x)dx for all polynomials p : R n — * R of degree at most 2. Then 



f(x)dx - T{f) 



^K/lloo / \\x-pfdx 

's 



foranyf£C 2 (S). 



Proof. Let g(x) 



_ IK/II 



6 S, By Lemma [5] the functions g + f and g — f are 



convex. Observe that the operator vo u S \ T fulfils the assumptions of Theorem [2] (for the 
purposes of this proof it is enough to have T defined on C 2 (S), the assertion of Theorem[2] 
holds also in this setting). Then 

vol(S) (g + f) (J!) < T(g + /) and vol(S) (g - /) (p) ^T(g-f). 

Using the first of the above inequalities, Lemma [6] and the second assertion of Lemma [5] 
we infer that 



1 



(g + f)(x)dx^vol(S)(g(p) + f(p)) + -\\d 2 (g + f)\\ oQ / ||x-p|| 2 ^ 
s z Js 



<T(«7 + /) + ||d 2 /||oc / l^-pll 2 ^ 



Since J„ g(x)dx = T(g), then 



(5) 



f{x)dx-T(f)^\\d 2 f\\ 



v I N x -p|| dx. 
s 



Repeating the above argument for g — f we arrive at 



f(x)dx - T(f) 



<||d 2 /IU [ \\x-p\\ 2 dx, 
Js 



which, together with ©, finishes the proof. 



□ 



To approximately compute multiple integrals on simplices using the cubature rules it is 
enough to have error bounds of such rules on the unit simplex Su i.e. on the simplex with 
vertices (0, 0, ... , 0), (1,0,..., 0), (0, 1, . . . , 0), . . . , (0, 0, . . . , 1). The integral over any 
simplex by the appropriate affine variable interchange can be computed as an integral over 
Si . That is why we record below the following 

Corollary 8. Let T be a positive linear functional defined on C 2 {S\) such that T{p) = 
f s p(x)dx for all polynomials p : R™ — » R of degree at most 2. Then 



f(x)dx - T(f) 



Si 



< 



(n + 2)l(n+ 1) 



foranyfeC 2 ^). 



Proof. It is well known that vol(5i ) = -^ . To prove the Corollary by Theorem [7] it is 
enough to check that 



\x — p\\ dx 



Si 



(n + 2)!(n+l) : 



where p 



n + 1 



n + 1 
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For x = (xi, . . . , x n ) € S\ we have 

-\ 2 



n 

\x-p\\ 2 = Y. 
t=i 



i 



i 



n r 

E 

z=i 



*?- 



1 



E 
j=i 



*?(*) 



(n+1) 2 
2 



n 



1 



m(x) 



i 



(n + 1) 2 



where 7T,(x) = x. L is the projection to the i-th axis. Bessenyei [ 1 , Lemma 1] computed 

1 



Repeating his argument we can easily compute 



1. 



, n. 



Trf(x)dx 



Then 



/ \\x-p\\ 2 dx = J2f 
JSi ,-_-, JS 



E 

2=1 



2 



(n + 2)!' 

2 



71+ 1 

2 



i« (a:) 



(n+1) 2 



dx 



1 1 

(n + 2)! ~ n+1 ' (n + 1)! + {n + l) 2 ' n! 

n 2 



(n + 2)! (n + 1)! (n + 1) 



(n + 2)!(n+l)' 



D 



The above Corollary is useful in the approximate integration since an important class 
of positive linear functionals which coincide with the integral for polynomials of degree at 
most 2 are operators of the form of a conical combination of some values of an integrand 
at appropriately chosen points. Below we give an example of such cubature operator. 

Corollary 9. Let T be the unit simplex in R 2 , 

nf) ■■= y a (m °) + m !) + /c 1 ' °)) ■+■ §/ (I' I) • 



Then 



f(x,y)dxdy -T(f) 



<^f\U 



foranyfeC 2 (T). 

Proof. A cubature operator T fulfils all the assumptions of Corollary [8] in particular, 
T(p) = JJ T p(x, y)dxdy for any polynomial p : M 2 — > R of degree at most 2. □ 
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